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Abstract 

In  two  dimensions,  the  trajectory  of  a  constant  veloc¬ 
ity  missile  which  intercepts  a  zero-acceleration  target  and 
which  minimizes  the  weighted  sum  of  the  final  time  aniT 
the  integral  of  the  missile  normal  acceleration  squared  is 
derived.  The  launch  angle  of  the  missile  is  arbitrarily  pre¬ 
scribed.  The  optimal  control  involves  elliptic  functions  and 
requires  the  numerical  solution  of  two  non-linear  algebraic 
equations  foi  its  calculation.  This  procedure  is  used  to  cal¬ 
culate  the  time-to-go  in  a  six-degree-of-freedom  simulation 
in  which  the  target  performs  accelerated  maneuvers.  The 
current  missile  velocity  is  projected  on  the  plane  of  the  line 
of  sight  vector  and  the  target  vector,  and  the  time  required 
to  perform  the  two-dimensional  intercept  is  calculated  and 
used  as  the  time-to-go  for  the  linear-quadratic  guidance  law 
of  the  missile.  Results  show  substantial  improvement  rel¬ 
ative  to  the  range-over-closing-speed  method  of  computing 
time-to-go. 

List  of  Symbols 


m 


A.B.C 

C 


Missile  normal  acceleration  (ft/sec^) 

Constants  defined  in  Eqs.  (21)  through  (23) 

End-point  function 

Hamiltonian 

Performance  index 

Modulus  of  elliptic  functions 

Magnitude  of  missile  acceleration  oscillation 

Time  (sec) 

Ratio  of  target  speed  to  missile  speed 
Velocity  (ft/sec) 

Weight  in  the  performance  index 
Inertial  coordinates 

Nondimensional  missile  normal  acceleration 
Angle  of  maximum  normal  acceleration 
Nondimensional  relative  abscissa 
Nondimensional  relative  ordinate 
Orientation  of  missile  velocity  vector  to  X-axis 
Variable  Lagrange  multiplier 
Nondimensional  time 
Orientation  of  target  velocity  to  X-axis 
TVansformed  variable  representing  8 
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Introduction 


A  guidance  law  of  current  interest  for  bank-to-turn  hom¬ 
ing  missiles  is  the  linear-quadratic  guidance  law  which  con¬ 
tains  proportional  navigation  as  a  particular  case  (see,  for 
example.  Ref.  1).  In  order  to  implement  this  guidance  law. 
an  algorithm  for  predicting  time-to-to  is  needed.  The  sim¬ 
plest  time-to-go  formula  is  range  divided  by  closing  speed 
and  is  valid  for  a  constant- velocity  missile  and  target  on  a 
collision  course.  This  formula  has  been  improved  in  Ref.  2 
by  accounting  for  the  missile  longitudinal  acceleration. 

Unfortunately,  the  linear-quadratic  guidance  law  tends 
to  drive  the  missile  and  the  target  into  a  homing  triangle 
in  which  range  and  closing  speed  become  unobservable.  In 
Ref.  1,  a  linear-quadratic  guidance  law  for  dual  control 
(intercept  and  estimation  enhancement)  has  been  proposed. 
This  guidance  rule  moves  the  missile  away  from  the  homing 
triangle  improving  estimation  but  making  the  time-to-go 
algorithm  invalid. 

The  purpose  of  this  study  is  to  develop  a  time-to-go 
algorithm  which  is  suitable  for  dual  control.  This  is  accom¬ 
plished  by  finding  the  normal  acceleration  history  which 
minimizes  the  time  for  a  constant-velocity  missile  to  inter¬ 
cept  a  nonaccelerating  target  and  using  this  time  as  the 
time-to-go.  To  account  for  a  limit  on  the  normal  accelera¬ 
tion,  the  performance  index  is  weighted  by  the  integral  of 
the  control  squared. 

As  a  first  step,  the  analysis  is  carried  out  in  two  di¬ 
mensions  with  the  hope  that  some  insight  in  the  three- 
dimensional  problem  will  be  achieved.  The  resulting  time- 
to-go  algorithm  is  tested  in  a  six-degree-of-freedom  simula¬ 
tion  by  projecting  the  current  missile  velocity  vector  onto 
the  plane  of  the  line  of  sight  and  the  target  velocity.  The 
optimal  intercept  time  is  computed  in  this  plane  and  used  as 
the  time-to-go  for  the  linear-quadratic  guidance  law.  Missile 
acceleration  is  included  by  using  an  average  missile  speed. 

Optimal  Intercept  Problem 

The  geometry  of  the  intercept  problem  is  shown  in  Fig.  1 . 
The  XY  coordinate  system  represents  an  inertial  frame,  and 
the  X  axis  is  along  the  line  of  sight  at  t  =  0.  The  target, 
located  at  Xt  =  Jfo  at  1  =  0,  is  moving  along  a  straight  line 
which  makes  an  angle  ^  with  respect  to  the  X  axis.  The 
constant-speed  missile  is  launched  at  an  angle  relative  to 
the  A'  axis,  and  the  velocity  direction  S(<)  is  changed  by 
controlling  the  normal  acceleration  a„{t). 


Ill  ternns  of  the  nondimensional  quantities 

_  Vm  ,  Xt-Xm  _  Yt-Yk, 

^  ^  Y  ^  '  —  Y  ’  ^  “  Y 

0  ^  o  ^  o 


Vg,'  V*,’ 

the  optimal  intercept  problem  is  stated  as  follows:  Find 
the  normal  acceleration  history  Q„(r)  which  minimizes  the 
performance  index 

.  ^  l-W  ri  ,  . 


subject  to  the  differentia]  constraints 

^  —  V  cos  <ti  —  cos  8  (3) 

jj  =  t)sin0  — sintf  (4) 

e  =  a„  (5) 

cind  the  prescribed  boundaiy  conditions 

=  0,  fa  =  1,  i;o  =  0,  tfo  =  <?os  (6) 

0  =  0,  =  0  (7) 

where  Oo,  is  the  specified  value  of 

The  performance  index  is  the  weighted  sum  of  the  final 
time  and  the  integral  of  the  normal  acceleration  squared. 
Without  the  integral,  the  solution  of  the  optimal  control 
problem  is  infinite  normal  acceleration  to  rotate  the  missile 
velocity  to  the  homing  triangle  direction  followed  by  zero 
normal  acceleration  until  intercept.  The  physical  bound 
which  exists  on  normal  acceleration  can  be  incorporated 
directly  as  a  control  variable  inequality  constraint.  This 
approach  makes  the  optimal  control  be  maximum  normal 
acceleration  followed  by  zero  normal  acceleration  and  is 
the  subject  of  a  companion  study  (Ref.  3).  Here,  the  nor¬ 
mal  acceleration  bound  is  included  indirectly  by  imposing 
a  penalty  on  the  performance  index  for  high  normal  accel¬ 
eration.  By  adjusting  the  weight  \V,  the  highest  normal 
acceleration  encountered  along  the  trajectory  can  be  con¬ 
trolled. 

Solution  of  the  Optimal  Control  Problem 

From  the  variational  Hamiltonian 
1  —  IF 

H  =  — - — o’ -t- A((i'cosd  -  cosfl) 


From  the  end-point  function 

G=W(Tf-  r.)  +  -t-  (13) 

it  is  seen  that  the  natural  boundary  conditions  are 

H,  =  -G,,  =  -W  (14) 

~  *'t 

X,j  =  G„=v,  (16) 

=  G,,=Q  (17) 

Eqs.  (12)  and(17)  imply  that 

a„/=0  (IS) 

Since  the  Hamiltonian  is  not  an  explicit  function  of  time, 
the  first  integral  H  =  const  exists.  When  combined  with 
Elq.  (14),  the  first  integral  becomes 

II  =  -W  (19) 

Hence,  if  the  definition  (S)  of  H  is  substituted  into  Eq. 
(19),  the  following  expression  for  the  optimal  control  can 
be  obtained 


On  =  ±\ 


-{A  +  BsiitO  +  CcosO) 


A  =  IF  •+ Ajiicos^ -f  A,usin  0  (21) 

B  =  -An  (22) 

C  =  -A('  (23) 

If  0,  were  free,  the  optimal  a,  would  be  o,  =  0,  and  the 
optimal  intercept  path  would  be  =  const.  =  9sl  or  a 
straight  line.  Hence,  if  0^  <  Osi,  the  plus  sign  in  Eq.  (20j 
holds,  and  if  Bp  >  Bsl,  the  minus  sign  is  valid. 

Eq.  (20)  can  also  be  obtained  by  integrating  Eq.  (11). 
This  is  accomplished  by  changing  the  independent  variable 
from  T  to  using  Eq.  (5),  using  Eq.  (12)  to  relate  a„  to  A#, 
integrating,  and  applying  the  final  condition  (17).  Finally, 
9/  is  eliminated  by  using  the  final  condition  (14). 

In  order  to  integrate  the  differential  equations  (3)  through 
(5),  it  is  useful  to  rewrite  Eq.  (20)  in  the  elliptic  function 
form  (Ref.  4) 


On  = 


(A  +  P){\-K^s\nU') 


-1-A,(usind  -  sind)  +  As(a„), 

it  is  seen  that  the  Euler-Lagrange  equations  are  given  by 

A(  =  -tU  =  0  (9) 

A,  =  -//,  =  0  (10) 

=  -//*  = -A(sintf-(- A,costf  (11) 

0  =  =  (1  -lF)a„-l-A,  (12) 

Eqs.  (9)  and  (10)  state  that  A(  and  A,,  must  be  constants. 
Before  solving  the  remaining  equations,  the  natural  bound¬ 
ary  conditions  are  derived. 


tan7=§.  P  =  s/B^  +  C\K^  =  ^  (25) 

If  T  is  chosen  to  be  the  dependent  variable,  Eq.  (5)  can 
be  integrated  as 

dO  f*i  dxi, 

r,p  =  Tf  -Tp=  /  — ^  =  2  — — -  (26) 

Jip  Qn(<l)  o„(v) 

Then,  substitution  of  Eq.  (24)  leads  to 


V  =  ±/<'\ 


where  F(rl>,  K)  is  the  elliptic  integral  of  the  first  kind. 
Similarly,  Eqs.  (3)  and  (4)  can  be  expressed  as 


1-Vo 


[v  ain 


(28) 

(29) 


Carrying  out  the  integration  and  applying  the  initial  condi¬ 
tions  (6)  leads  to 


-^1E(*.a')-e(*,.a)1  (30) 


and 


.7=±K  «n  ^  )  .in  -»)(#■  (*,K)-F(*„A')1 

-i^(E(*.A)-£(*,.A)l  (31) 

+  A'*  iin^  .in^ 

where  £(^,  K)  is  the  elliptic  integral  of  the  second  kind. 

The  last  step  is  to  satisfy  the  final  conditions  (7).  This 
amounts  to  finding  the  values  of  X(  and  A,,  which  satisfy 
Eqs.  (30)  and  (31)  evaluated  at  t/>y  where  =  q/  =  0.  Be¬ 
cause  these  equations  are  nonlinear,  the  solution  process  is 
iterative.  Rather  than  solving  for  the  Lagrange  multipliers, 
it  is  more  convenient  to  work  in  terms  of  and  0/.  The 
relationships  between  the  t!>’s  and  the  A’s  are  obtained  from 
Eqs.  (24)  and  are  given  by 


A, 


-t4^cos(tfo -2t^p) 

2sin’t(’/  -  I  +  vcos{(i)  -0o  +  20„)  '  ' 

-lk'sin(g, -20,) 

2sin^  0/  —  1  -f  t/cos(0  -0„  +  20,)  '  ' 


Numerical  Results 


weight  has  been  chosen  to  be  IV  =  0.4.  This  value  gives 
reasonable  trajectories  for  all  parameter  values. 

The  figures  presented  here  have  been  generated  mainly 
to  verify  the  computational  procedure.  Figs.  2  through  4 
contain  0,  and  0/.  Next,  the  normal  acceleration  histories 
are  shown  in  Figs.  5  through  7  for  the  case  where  v  =  0.5 
or  Vm  =  2Vt-  Then,  Figs.  8  through  10  present  the  flight 
time,  and  finally,  the  trajectories  are  presented  in  Figs.  1 1 
through  13.  The  results  presented  in  all  the  figures  seem 
reasonable  in  that  they  follow  expected  trends. 

Time-To-Go  Prediction  for  an  Accelerating  Missile 

It  is  desired  to  use  the  previous  results  in  a  six-degree- 
of-freedom  homing  missile  simulation  in  which  the  missile 
velocity  is  not  constant  and  the  target  acceleration  is  not 
zero.  First,  the  missile  velocity  vector  in  the  simulation 
is  projected  onto  the  plane  of  the  line-of-sight  vector  and 
the  target  velocity  vector,  and  the  time-to-go  is  calculated 
in  this  plane.  Second,  a  nominal  tangential  acceleration 
history  is  assumed  for  the  missile  and  a  constant  average 
missile  velocity  is  determined  for  a  given  initial  time  and  an 
assumed  final  time.  This  velocity  is  then  used  to  compute 
the  time  for  intercept  as  defined  in  Section  3.  Since  the 
computed  final  time  and  the  assumed  final  time  are  not 
equal,  an  iterative  process  is  used  to  compute  the  final  time 
and,  hence,  the  current  time-to-go. 

Some  results  are  presented  in  Table  1  for  a  target  which 
performs  two  maximum  acceleration  maneuvers;  one  when 
the  missile  is  6,000  feet  away  and  one  when  the  missile  is 
1.0  sec  away.  Miss  distances  are  shown  for  two  ofT-boresight 
angles,  two  launch  ranges,  and  seven  aspect  angles.  In  Case 
I,  time-to-go  is  calculated  as  range  over  closing  speed,  and  in 
Case  II,  it  is  calculated  using  the  minimum-time  solution.  In 
general,  the  new  time-to-go  algorithm  yields  a  smaller  miss 
distance  than  range  over  closing  speed.  In  several  cases  it 
gives  a  hit  where  the  latter  gives  a  miss  (more  than  ten 
feet). 

Discussion  and  Conclusions 


Given  values  for  W,  v,  0,  and  the  values  of  0,  and 
0/  which  satisfy  Eqs.  (30)  and  (31)  evaluated  at  the  final 
point  are  computed  using  Newton’s  method  with  analytical 
derivatives.  Given  0,  and  0/,  A;  and  A,  follow  from  Eiqs. 
(32)  and  (33).  Next,  Eqs.  (21)  through  (23)  give  A,  B,  and 
C,  and  Eqs.  (25)  give  7,  P,  and  K.  Finally,  the  values  of 
and  T)j  can  be  computed  from  Eqs.  (30)  and  (31).  Improved 
values  for  v,  and  0/  (values  which  make  and  17/  closer  to 
zero)  are  provided  by  Newton’s  method. 

Numerical  results  have  been  generated  for  the  following 
parameter  values; 

u  =  0,  0.2,  0.4,  0.6,  0.8 
0  =  0",  45',  90' 

-70'  <  tf,  <  70' 

The  limits  on  fl,  ate  dictated  by  the  field  of  view  of  the 
seeker.  Results  are  available  for  0  up  to  180'in  Ref.  4, 
but  they  have  not  been  presented  here  because  of  space 
limitations.  After  some  experimentation,  the  value  of  the 


A  time-to-go  algorithm  has  been  developed  which  is 
valid  for  missile- target  geometries  which  differ  greatly  from 
the  intercept  triangle.  The  algorithm  is  based  on  the  analyt¬ 
ical  solution  of  a  weighted  minimum-time  problem.  While 
the  algorithm  has  been  developed  for  a  constant  velocity 
missile  in  two  dimensions,  it  can  be  applied  to  an  acceler¬ 
ating  missile  in  three  directions  by  using  an  average  missile 
velocity  and  by  computing  the  time-to-go  in  the  plane  of 
the  line-of-sight  vector  and  the  target  velocity  vector.  A 
six-degree-of-freedom  simulation  shows  that  the  proposed 
time-to-go  algorithm  improves  miss-distance  performance 
relative  to  computing  time-to-go  as  range  divided  by  clos¬ 
ing  speed. 
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Table  1:  Simulation  Results 


Off. 

Ooresight 

Angle 

(deg) 

Launch 

Range 

(ft) 

Aspect* 

Angle 

(deg) 

Case  I 
Miss 
Distance 
(ft) 

Case  11 
Miss 
Distance 

(ft) 

0 

7000 

0 

0.156 

0.337 

0 

7000 

30 

0  694 

0.243 

0 

7000 

60 

1.203 

0.284 

0 

7000 

90 

1.322 

0.469 

0 

7000 

120 

2.561 

0.265 

0 

7000 

150 

0.739 

0.077 

0 

7000 

180 

0.351 

0.010 

0 

3000 

0 

0.983 

0.097 

0 

3000 

30 

0  703 

0.069 

0 

3000 

60 

1.382 

0.063 

0 

3000 

90 

2.512 

0  138 

0 

3000 

120 

11.85 

2.950 

0 

3000 

150 

4.166 

2911 

0 

3000 

180 

0.605 

0.253 

•to 

7000 

0 

301.7 

1251. 
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Fig.  1  ,  Intercept  Geometry 
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Fig.  2a,  i/>„  for  ^  =  0  deg 


PHI=  0  DEG. 


Fig.  2b,  1/7  fof  0  =  0  deg 
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Fig.  1 1  ,  Trajectories  for  p  =  0  deg 


PHI=  45  DEG. 
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